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Abstract 

It has been known for 30 years that 't Hooft-Polyakov monopoles of charge Q greater 
than one cannot be spherically symmetric. 5 years ago, Bolognesi conjectured that, at 
some point in their moduli space, BPS monopoles can become approximately spherically 
symmetric in the high Q limit. In this note we determine the sense in which this conjecture 
is correct. We consider an SU (2) gauge theory with an adjoint scalar field, and numerically 
find configurations with Q units of magnetic charge and a mass which is roughly linear 
in Q, for example in the case Q = 81 we present a configuration whose energy exceeds 
the BPS bound by about 54 percent. These approximate solutions are constructed by 
gluing together Q cones, each of which contains a single unit of magnetic charge. In each 
cone, the energy is largest in the core, and so a constant energy density surface contains 
Q peaks and thus resembles a sea urchin. We comment on some relations between a 
non-BPS cousin of these solutions and the dark matter halos of dwarf spherical galaxies. 
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1 Introduction 



Within the large moduli space of solutions of BPS monopoles [UEIEHI] of charge Q, it 
is plausible that there exists a high Q limiting sequence of monopoles that are spherically 
symmetric up to 1/Q corrections. However even the most qualitative features of such solu- 
tions are, as yet, unknown. In this note we will solve an easier problem, we will explicitly 
construct configurations in an SU(2) gauge theory with an adjoint scalar whose energy 
slightly exceeds the BPS bound. These configurations do not provide time-independent 
solutions of the equations of motion^, but the near saturation of the BPS bound supports 
a conjecture that certain qualitative features of these solutions will be shared by some true 
solutions. In particular, our approximate solutions become spherically symmetric in the 
high Q limit, and so the corresponding true solutions are also asymptotically spherically 
symmetric. 

We will construct these approximate solutions by decomposing space into Q identical 
cones which extend from the origin, each of which asymptotically contains a unit of 
magnetic flux and is axially symmetric in a sense which will be made precise below. Of 
course, the classification of regular polyhedra implies that no decomposition of space into 
identical cones exists for general Q. Therefore the cones may not fill all of space. We will 
provide configurations in which the space between the cones consists of a vanishing gauge 
field and a continuous Higgs field, which at high Q provide a contribution to the energy 
which decreases with Q, and so is subdominant with respect to the total energy of the 
configurations, which is proportional to Q. 

Our problem is then reduced to the following steps. First, in Sec. |2]we will choose 
an Ansatz and boundary conditions for the configuration in a single cone. The finiteness 
of the energy fixes all of the functions in our Ansatz except for two angular functions 
which must be chosen. Different choices yield different energies, none of which satisfy 
the BPS bound but all of which satisfy the bound up to a factor which is Q-independent 
at large Q. This Q-independence, which is critical for the qualitative agreement with a 
BPS configuration, is achieved if we impose that the radial dependence of the solution 
satisfies an ODE in a single variable. Next in Sec. |3]we will explicitly solve this equation 
in the large and small radius limits, and discuss how these solutions grow together. In 
this way we are able to establish the qualitative profiles of these solutions, confirming the 

1 Thcy do provide initial conditions for oscillatory time-dependent monopole solutions. 
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monopole bag description conjectured in Ref. j5]. We will then present numerical solutions 
of the ODE in Sec. HI Given the large moduli space of solutions, one may wonder why 
approximately spherically symmetric monopoles are interesting, especially given the fact 
that exact solutions which are spherically symmetric up to 1/Q corrections are unknown. 
In Sec. |5]we describe the kinds of field theories in which we expect such monopoles to be 
the only monopoles which survive a non-BPS deformation, and provide a very speculative 
application of such theories. 

Since Bolognesi's groundbreaking proposal [5j, the field of monopole bags has been 
steadily advancing with the understanding of various limiting behaviors [6] and of large Q 
limits of Nahm transforms of bag solutions [7j. Yesterday a paper by Manton appeared on 
the arXiv with significant overlap with our results [8]. He found exactly BPS Bolognesi bag 
solutions with an approximate spherical symmetry of the type considered here. This begs 
the question: how we can justify presenting an approximate solution after exact solutions 
have appeared? Our justification is that the physical application that we have in mind in 
Sec. |5] is in a theory in which we expect the radial dependence to be very distorted, but 
under a radial deformation the cone structure of our configurations is preserved. 

2 Ansatz and equations of motion 



As described in the introduction, charge Q magnetic monopole configurations will be 
constructed from Q identical cones extending from the origin together with an extrapo- 
lation in the region between the cones which yields a contribution to the energy which is 
sub dominant at large Q. The crucial step in the construction is therefore to provide the 
configuration in a single cone. The first step is to determine the size of the cone and to 
fit it with coordinates. 

For simplicity let Q = n 2 be a perfect square. This limitation will not be relevant for 
our analysis of the large Q regime. We will consider a cone whose core extends along the 
positive z axis, and will let ip be the azimuthal angle 



2.1 The cone 



tan^ = — . 

x 



(2.1) 



If we define the radius in cylindrical coordinates to be 



p = y / x 2 + y 2 , 



(2.2) 
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then the cone will be defined by the condition 

P < o- , (2.3) 

n 

where the cone is small enough that n 2 non-overlapping copies fit inside of M 3 . The 
constant a parametrizes the size of the cone. 



The base of each cone at a radius r = a/x 2 + y 2 + z 2 is then a circle of radius p = 
r I + 1 ~ or jn and so the area of the base of the cone is about na 2 r 2 /n 2 = -na 2 r 2 /Q 
which is cr 2 /4<5 times the total area at that radius, implying that o < 2 in order for the 
area of the cones to be less or equal to the area of the spatial sphere. This is still not 
sufficient for the Q cones to fit inside of M 3 . Here we are interested in the large Q limit 
and hence we can consider the size of the base of the cone to be much much smaller than 
the radius of curvature of the sphere at r. Then the limit, which is an upper bound on 
the size of the cone, is given by the size of the circle which can fit inside of a hexago 
More precisely, in the large Q limit the upper bound on a is 



<W = 2 J-^- = 2 A /^-l-9. (2.4) 

^hexagon y Zi \ O 

Clearly such a configuration cannot hope to saturate the BPS bound, instead a qual- 
itative agreement with a true solution leads to the requirement that it exceeds the BPS 
bound by a factor (of order unity) which tends to a constant in the high Q limit. 



2.2 The Ansatz 



We will consider an £77(2) gauge theory with a massless adjoint scalar $. Furthermore, 
we will make the crude approximation that the configuration of the Higgs and gauge field 
factorizes into z and p/z-dependent functions, yielding the following axially-symmetric 
Ansatz 



$ = h(z) F(rj) (c^ + st 2 ) -F 2 (r])t 
for the Higgs field and 

a(z) 



(2.5) 



A 1 
A 3 



z 
a(z) 

z 
a(z) 
nz 



(cs [J{rj) - G( V )} t 1 + [c 2 G(r]) + s 2 J{r 1 )] t 2 - sH( V ) t 3 ) , 
(- [c 2 J( V ) + s 2 G( V )] t 1 - cs [J( V ) - G{r])} t 2 + cH( V ) t 3 ) 
I^ist'-ct 2 ) , 



(2.6) 



2 This hexagonal lattice approximation is exact for monopoles in AdS [91 110). 



3 



for the 577(2) gauge field Ai where we have defined the variables 

TIP 

77 = — e [0, a] , e = sign(F'(?7)) , c = cos-0, s^sin^, (2.7) 

and the functions F, G, 77, 7 and J of i] and also ft, and a of 2. 

Such solutions are invariant under a rigid axial symmetry which simultaneously rotates 
the vectors 

cos0 — sin0\ fx t 1 
sin0 cos0 / \y t 2 



1 X 








\y 





(2.8) 



and the potential 

Ai \ I cos sin \ / Ai 
Ai) \ — sin cos (j) j \A 2 

by an arbitrary angle 0. 



(2.9) 



2.3 Finite energy conditions 
Topological conditions 

We define the charge of the magnetic monopole to be equal to the number of cones Q = n 2 
in which the Abelian magnetic field Tr(FQ) asymptotically pointing in the outward radial 
direction integrates, over the base of each cone, to a single Dirac quantum 4-7T. The 
configuration then will have a finite energy only if, for all sufficiently high z, the value of 
$ on the cone's base is in the correct topological sector. More precisely, we will demand 
that, at large z, $ be constant on the boundary p = az/n of the cone. Therefore the value 
of $ on the base of the cone defines a map from a 2-sphere S 2 , which is the (z — 00, rj < a) 
two-disc D 2 with its boundary identified with a point, to the space S 2 of values of $ of 
constant norm |$| = v. The topological condition is then that this map S 2 — > S 2 be of 
degree one. 

We will satisfy this condition by imposing that F be a continuous function such that 

0<F(?7)<1, F(0) = 0, F(cr) = 0, (2.10) 

and F(i]q) = 1 for some intermediate value i]q of its argument. We will assume that 
the function F has only one maximum on the interval [0, a]. One function obeying this 
condition is 

F(rj) =sin(^) . (2.11) 
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We will see shortly that the analysis of the region between the cones is simplified if one 
requires the derivative of F to be zero at 77 = a which can be achieved via the deformation 



where k 3> 1 is a large integer. With this deformation, F needs to be rescaled so that its 
maximal value is equal to one. It is difficult to calculate the resulting normalization of F 
analytically, but for k 3> 1 the normalization constant is arbitrarily close to unity. 

High z asymptotics 

The topological condition on the Higgs field is necessary but not sufficient for the finite 
energy of the configuration. The total energy is the sum of the integrals of the gauge field 
strength energy density \Fij\ 2 and the kinetic energy density of the Higgs field |-Dj$| 2 . So 
long as the local energy density is finite, a divergence may only arise from the noncompact 
region z — > 00, where we will demand that 



and so their derivatives tend to zero. The former condition is an arbitrary choice, however 
in the non-BPS generalizations discussed in Sec.[5]it minimizes the potential for the Higgs 
field. As will be clear from the analysis that follows, any different choice of limiting value 
for a will lead to a rescaled value of the finite energy conditions on the functions G, 
H, I and J. The combination of the rescaling of a with that of the functions leads to 
precisely the same form of the connection, and so this reflects a simple redundancy in the 
parametrization of our Ansatz. 

Higgs field kinetic energy 

We will now impose that the total energy in each cone is finite. There are two contributions 
to the energy, one from the magnetic field F^F^ and another from the kinetic term of the 
Higgs field |_Dj$| 2 . As both contributions are positive definite, we must demand that they 
are finite separately. We now begin by imposing the finiteness of the Higgs field kinetic 
energy. 

The finiteness of the Higgs field kinetic energy places non-trivial conditions on the 
functions in our Ansatz. At high z the Higgs field tends to a constant limit, and so its 
covariant derivative is at most of order 0(1/ z) and the energy density |-Dj<3>| 2 is therefore at 
most of order 0(1/ z 2 ). The finiteness of the energy of the configuration is then equivalent 




(2.12) 



h(oo) = gv 



a(oo) = n , 



(2.13) 
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to the vanishing of the order 0(l/z 2 ) and 0(1/ z 3 ) terms of |-Dj$| 2 for each spatial index i 
and each gauge direction. As there are three spatial directions and three gauge directions, 
this consists of 9 directions, of which we will see 3 are independent. 

The axial symmetry ( 12. 8 p and (I2.9p implies that the energy density is axially symmet- 
ric. Therefore it suffices to calculate the energy at tp — 0, corresponding to x = p and 
y = 0. This is merely to simplify the expressions in the following. The three contributions 
to the kinetic energy are then the squares of the quantities -Dj$ = <9j$ + i[Ai, $], 

D x $\ y=0 = ^(n\F'\-aGVT^^--^L=t^ , 

Dy^o = v(^ a ( eJ ^y^ +if )) t2 ' (2 - 14) 

The terms of order 0(1/ z 2 ) and 0(1/ z 3 ) in the kinetic energy will only arise if there are 
terms of order 0(1/ z) in _Dj$. Therefore if we assume that h — gv and a — n go to zero 
at least as quickly as 1/z, then we may drop the first term in the expression for -D z $ and 
approximate a by n everywhere in ( I2.14p without affecting the divergent terms. 

The 0(1/ z) terms in each of the three _Dj$ can then be seen to be proportional to a 
combination of the functions F, G, H, I and J. By setting these combinations to zero, 
we eliminate the divergence. These three conditions can be solved to yield, for example, 
G, H and / as functions of F and J, we find respectively 

\F'\ 1 Vl - F 2 n\F'\ 

G= ' 1 , H = --e. r ' , 1= J' 1 1 . (2.15) 

Recall that F(a) = therefore H is only finite if J(cx) = as well, and in fact the ratio 
of the two needs to tend to a constant as z tends to a. F(0) = as well, and so the 
finiteness of H(0) requires that the divergences in both terms cancel. Here e = 1 and so 

lim = lim i + 0(1) . (2.16) 

Substituting these relations into ( I2.14p the J dependence cancels and so we can find 

exact expressions for the covariant derivatives as functionals of F alone 

hF' , / , eF , 
D x $L- = (n-a) [ t 1 t 3 

hF 

D v <f>\ v=0 = —(n-a)t 2 , (2.17) 

y y ZT] 



D z $\ y=0 = h'iFt 1 + eVl-F 2 t 3 ) - ^D x $ 



\y=o 
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Notice that finite energy densities require 

i h , F , 

lim — < oo , lim — < oo , (2-18) 

z-s>0 Z V^O Tj 

and imply in particular 

h(0) = a(0) = , (2.19) 

as they must be in order for $ and to be well defined at the origin. A finite energy 
density at the maximum of F also requires that F' vanish as quickly as yl — F 2 . 

The total Higgs kinetic energy density is found by adding the squares of the compo- 
nents in Eq. d2~T7jl 

^ = ^E^»i' = ^(V)' + ^(»-<.)»(5 + S( 1 + 8) • (2 ' 20) 

As the covariant derivatives were independent of J, so is the total kinetic energy density. 
Magnetic field kinetic energy 

As was the case for the Higgs kinetic term, the axial symmetry of the energy means that, 
for the purposes of calculating energy, it suffices to consider the gauge field strength 

F i:i = diAj - djAi + i[A h Aj] , (2.21) 

at y = 0. Again this contribution to the total energy will be finite if the energy density 
is at most of order 0(1/ z 4 ), which means that the field strength itself must be at most 
of order 0(1/ z 2 ). Therefore it will again suffice to impose that terms of order 0(1/ z) 
vanish, and that the coefficient of the order 0(1/ z 2 ) terms is finite. 

These conditions are only nontrivial in the case of F\2, which is 



'Vi\y=0 - ^2 



n 



./' • ^— ^ ) - (\GH 

V 



t' + 

A divergence at the tip of the cone can be avoided if 



n ( H' + — ) - aGJ 



V 



t 3 ) . (2.22) 



lim — < oo . (2.23) 

The expression f)2.22p is then, at any constant 77, of order 0(1/ z 2 ) and so the energy is 
only divergent if the field strength itself diverges. To avoid such a divergence as 77 tends 
to we will impose that the two quotients by 77 are finite 



hm < 00 , hm — = hm — < 00 , (2.24) 

ry— >0 fj V— >® Tj »?->0 \Tj TjF 
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where we have used the expression for H in Eq. ( I2.15p . From here we learn that G(0) = 
J(0), H(0) = and also this yields a refinement of the boundary condition ( I2.16P 

lim — = lim - + 0(n) . (2.25) 

If we now expand F as 

F ~ Cj-rf + O (r] j+1 ) , (2.26) 

where j = 1,2,3,... and j represents the first non-zero term in the expansion then 
Eq. (12T25]) implies 

J ~ Cjrf' 1 + O (f] j+1 ) . (2.27) 

This is sufficient for rendering H/n finite at small rj. Finally, G ~ jcjrf' 1 + • • ■ together 
with Eq. (I2.27P implies that the first constraint of Eq. (I2.24p is trivially satisfied for j > 1 
while it is also non-trivially satisfied for j = 1 due to the matching of the coefficients of 
G and J. 

The other components of the field strength can be evaluated easily. Using the fact 
that / = rjG from Eq. ( 12. 15ft one finds the y component of the magnetic field 

Fi2,\ y= Q = ——Gt 2 , (2.28) 

which yields a finite energy contribution at large z since a 1 is at most of order 0(l/z 2 ). 
A divergence is avoided at small z if in addition to (I2.23P one imposes 

a' 

lim — < oo . (2.29) 

z->0 Z 

In fact this condition follows from (I2.23|) if a is differentiate at r\ — 0. 
The final component of the magnetic field is 





~ar] 1 


\y=o - 


_z 2 \ 



ar) ( G — J aGH j\ a' J 

•q n 



t 1 



an ( H aJG rr \ a'H 
-4 — + H') + 



z 2 \ rj n J z 



t 3 . 



(2.30) 

Again the finiteness of this contribution to the energy is guaranteed by the condition 
Eq. d223}. 



2.4 What lies outside of the cones 

To choose a consistent configuration outside of the cones, one must first determine the 
configuration on the boundaries of the cones. We have already imposed the boundary 
condition F(cr) = and we have seen that the finiteness of H(a) implies that J(cr) = 
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as well. In general the shape of the region between the cones is quite complicated. 
However, we are only interested in an approximately BPS configuration, which we define 
as a series of configurations at various values of Q such that at large Q the energy is 
asymptotically proportional to Q. We will see in this subsection that it is therefore 
sufficient to consider a configuration in which the gauge field vanishes in the region between 
the cones, considerably simplifying our analysis. 

The vanishing of the gauge field in the region between the cones means that it also 
vanishes on the boundaries of the cones 

G(a) = H(a) = 1(a) = . (2.31) 

According to Eq. (I2.15p . for G(a) and 1(a) to vanish it is sufficient to fix 

F'(a) = , (2.32) 

which is the reason for the modification proposed in Eq. (12.121) . The vanishing of H(a) is 
only slightly more complicated. By (I2.15p . now with e = — 1, 



_ 1 = llm . (2 . 33 ) 

Since F'(a) vanishes, so must J' (a) and so taking yet another derivative of the numerator 
and denominator 

J'(a) = , J" (a) = . (2.34) 

a 

With these choices the gauge field vanishes on the boundary and the Higgs field is 
equal to 

$ = -h(z) t 3 = -h (r/^fi + a 2 /n 2 ^j t 3 . (2.35) 

These fields may then be easily extended to the region outside of the cones, simply by 
asserting that Aj always vanishes and $ always obeys (12. 35 p . This choice certainly does 
not minimize the energy, but as we will now argue, it yields a contribution to the energy 
which at large Q is increasingly sub dominant. 

Notice that the energy density between the cones is simply |<9j$| 2 which depends only 
upon r. Thus it is equal to a Q-dependent constant multiplied by the area of a sphere of 
radius 1 which is not within a cone. This area tends to a constant at large Q. In fact, at 
large Q the maximal number of cones that can be packed into a finite volume approaches 
the maximum packing of their circular cross-sections on a plane, which is given by the 
ratio of the area of a unit hexagon to a unit circle 2v / 3/tt. The area at unit radius outside 
of a cone is then 8v^3 — Air and hence, integrating over the radius 

(8^-471-) f . , i2 (8-V/3-47T) f 9/ „ o 
£out = ^-j 2 >- J drr 2 Tr \d^\ 2 = '- J drr 2 (h'(r)f , (2.36) 
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yields the energy outside of the cones. 

Therefore at large Q the geometric factor asymptotes to a Q- independent constant. 
Also at large Q the a/n in the argument in Eq. (I2.35P becomes negligible. Therefore 
the only source of Q-dependence is in h itself, which interpolates between and the Q- 
independent constant gv . We will argue in Sec. [3] that this interpolation occurs over a 
region of size of order 0(y/Q). Therefore one expects that the derivatives will decrease, 
and in particular the derivative squared energy density will scale as 1/Q. However the 
range of integration scales as y/Q, and so one expects the total contribution to the energy 
between the cones to decease as 1/a/Q, becoming ever subdominant as compared with 
the BPS energy which is proportional to Q itself. 

2.5 The total energy and BPS equations 

At an arbitrary point in space, the covariant derivatives of $ are 

Dx$ = (c 2 X 1 + s 2 X 2 ) t 1 + cs (X 1 - X 2 ) t 2 - cX 3 t 3 , (2.37) 
D 2 $ = cs {X x - X 2 ) t 1 + (s 2 X 1 + c 2 X 2 ) t 2 - sX 3 t 3 , (2.38) 
D 3 $ = [Z x - ^Ax) (ct 1 + st 2 ) + (Z 2 + ^A 3 ) t 3 , (2.39) 

where we have defined the following functions 

hF' hF , hFG , , « 

Ai = (n - a) , X 2 = — (n - a) , A 3 = (n-a) , 2.40 

Z ZTj Z 



Z x = h'F , Z 2 = h'eVl - F 2 , (2.41) 

and the field strength components are 

F 12 = -Ut (ct 1 + st 2 ) + -U 2 t 3 , (2.42) 

7] 7] 

F 23 = [c 2 (WJ + U x ) + s 2 WG] t 1 + cs (W(J -G) + U x ) t 2 + c (U 2 - WE) t 3 , (2.43) 

F 31 = cs (W(J -G) + C/i) t 1 + [c 2 WG + s 2 (WJ + U x )] t 2 + s (U 2 - WH) t 3 , (2.44) 



where we have defined 



TT ~ a 



G-J-rj(j'+ -Gh)} , U 2 = % \h + r)H' - %GJ 
\ n J J z 2 i n 



W 



a' 
(2.45) 
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The BPS equations read 



n n 
Z x - ±X X - -U x = , 
n 77 



77 n 
Z 2 + ±X 3 - -U 2 = 
n 77 



(2.46) 



= W J + C/i , X 2 = WG , X 3 = W# - C/ 2 . 
while the energy density for BPS-saturated configurations reads 
9 2 n conc , BPS = e^Tr (F jk $) = e^ k Ti (F^D^) 



n 



V 



{U X Z X + U 2 Z 2 ) + W (JX X + GX 2 + HX 3 



(2.47) 



The only function of rj which is not given in terms of the function F is J. 

We will in the following formally integrate the above boundary term, obtaining the 
correct magnetic charge. This procedure imposes some criteria on J. Using the relations 
( I2.15P we obtain the topological contribution to the energy density 



g 2 n c 



,BPS 



ah' 

Z 2 Tj 



. eFF' nriJF' nU + riJ') 
(n — a) — + 



F 2 



a'h n . 
+ —2{n-a) 

Z Z T) 



The total energy of the cone is given by 

2?r 
n 



E r 



dz I drj z 2 rj T-L c 



eFF' 
(2.48) 



(2.49) 



■'cone n 

' 2 J JO 

which tells us that it will be convenient to have the energy density of the form 
After this integration, the BPS part of the energy depends only on the boundary data 
and so by Stokes' theorem is of the form -^^f^-. Rewriting Eq. (I2.48P as discussed 



g 2 U 



cone,BPS 



+ ^2(n-«)^f-ev / T^F2 

^77 V 



(2.50) 



we see from the first and third terms that —t\/\ — F 2 \q = 2. Hence, in order to obtain a 
total derivative in z, we need to impose 



rjj 
"T 7 



(2.51) 



This can be done consistently with the boundary condition (I2.34p in many ways, among 
which we will choose 

F 

(2.52) 



J = -eVT^F 2 . 



V 
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Carrying out first the rj integration, we obtain 

f 1 r 2ah ' , \ 2a'h , v 2 r . ,, . , 

/ d77 77H(^e ) BPS = -^(2n-a) + -^2(n-a) = -^d z [(2n-a)toj , (2.53) 

which can readily be integrated 

4:7r Attv 

-E CO ne,BPS = -5— t [(2n - a)/^]^ = . (2.54) 

g z n z g 

This is exactly one Dirac quantum contained in a single cone. 

We will however see that we can only approximately satisfy the above BPS equations 
and hence we will need to calculate the total energy density 



H -J- 

'T-conc n n 

2g A 



X 2 + X 2 + X\ + [z x - ^X^j + (Z 2 + ^X 3 ) (2.55) 

+ + (Uf + U 2 2 ) + 2W (JUt - HU 2 ) + W^ 2 (G 2 + H 2 + J 2 ) , 

for our approximate solutions. The first five terms are the kinetic energy contributions 
given in Eq. fT2T2"U]) . 

The total energy of the monopole is then 

-^monopole Q^cone -^out • (2.56) 

Ideally the BPS equations could all be satisfied, in which case -E CO ne = -E C one,BPS, however, 
as we will see this does not turn out to be the case for our configuration. Hence the 
energy of the cone is given by Eq. (I2.49P with "H con e given by Eq. f)2.55p . This energy will 
necessarily exceed the BPS bound (I2.54p . 



2.6 Radial profile 

The conditions in Sees. 12.31 and 12.41 guarantee that the energy of these configurations 
is finite. However, a qualitative agreement with true BPS solutions is only possible if 
the total energy is, at large Q, proportional to Q. Such a scaling constrains the radial 
dependence of the Higgs field and gauge fields. There are many different choices of 
radial dependence which yield the correct scaling^. We will choose one of the simplest, 
we will impose that the Bogomol'nyi equations be satisfied near the cores of the cones. 
More precisely, we will expand the solution as a power series in 7]/n and will apply the 

3 Indeed two such choices were described yesterday in Ref. [S]. 
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Bogomol'nyi equations at leading order in the expansion. This will determine the radial 
profile functions. 

If we expand the function F as follows 

F(rj) =aq- hr] 3 + C(r/ 5 ) , (2.57) 
6 

and choose 



J = — eVl - F 2 , (2.58) 
rj 

as was suggested in Eq. (12.521) . we obtain at leading order in rj/n the following ODEs 

Hz) = -^L- , (2.59) 
n — a{z) 

h'(z) = ?-^(2n - a(z)) , (2.60) 
which can be combined into a single equation for a 

a"(z) ( a'(z) \ 2 a 2 a(z) . . /n ^\ 

~T~\ I P(2n-a z =0. (2.61) 

n — a{z) \n — a{z) J z A 

With the radial differential equation at hand we are now ready to write down all of the 
angular profile functions in the next subsection. 



2.7 Choosing an angular profile function 

Summarizing the construction of this section, given the functions F and J one may 
determine all of the angular functions in the Ansatz (12.5112. 6p . The radial functions on 
the other hand are determined by the Bogomol'nyi equations (12.591) and (I2.60p at leading 
order in rj/n. 

For instance, if F is given by Eq. (I2.12p and J by Eq. (I2.52p . then the other angular 
functions are given by 

— cos (22) (1 - 21) - ^ sin (22) 

G(v) * , == 

^/l-sin 2 (f) (1-21) 2 

whose rj/n expansion contains (12.611) with 

a = - . (2.63) 



F 2 (v) 



(2.62) 
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As we have mentioned, we will not be able to satisfy the BPS equations everywhere 
(or equivalently to all orders in the rj/n expansion) due to the fact that true BPS solutions 
do not exactly factorize. Hence, in order to measure the excess energy with respect to a 
true BPS configuration we need to calculate the total energy density as follows 



9 2 n conc = ~ti(zf + (^V(rj, a) + Q(rj, <r)j 



h(zf (n - a{z)f 



(2.64) 



l 



a 1 - 



;K(r},a) + -zSfaa 



rr 



a{zf (2n - a(z)f 



which should be integrated over the volume of the cone as in Eq. (12.491) while the total 
energy is given by Eq. (I2.56p . All the above functions V, Q, TZ and S are defined such 
that their integral f Q drj rjX = const is independent of cr, where X = V, Q, TZ, S. 

By evaluating the integrals over the angular functions numerically we can write the 
total energy for the cone 



2tt 

J^conc — o ri 



. a 



-a z h'(zY+ 1.24— + 3.11 
4 V n 2 



(2.65) 



3,1^ + 



6.02 1.24 



rr 



a(z) 2 (2n - a(z)f 



One learns from this expression that there is a competition in energy which is determined 
by the value of a. Only the first and the last term are dependent on a at large n and 
the first wants a to be small while the last one prefers a large value of a . It so happens 
that the last term wins the competition due to the fact that h varies only over a fairly 
small fraction of the integration range while the last term is roughly ^2 which is large 
for large n. Hence a minimization of energy leads to a value of o which is as large as 
possible, i.e. a ~ 1.9 as found in Eq. f)2.4p . We will check this statement numerically in 
Sec.H 



3 Asymptotics and qualitative features 

As we have described, the monopole is characterized by the configuration in a single 
cone. This consists of two parts, an angular function F which needs to be chosen as well 
as radial functions a, h. The factorizing Ansatz (I2.5f2.6p implies that the radial functions 
are independent of the choice of F. In particular the function a may be determined using 
the second order ODE (12.611) and then h may be determined from a using Eq. (I2.59p . 
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3.1 Small z asymptotic behavior 

At small z, near the center of the monopole, the gauge and Higgs fields approach zero. 
These two conditions imply the boundary conditions 

a(0) = a'(0) = . (3.1) 

In the small z region, a, a' <C n. Therefore (I2.6ip may be approximated by 

a»(z) - = • (3-2) 

z 1 

This is a homogeneous, linear, ODE with the solution 



a(z) = k lZ ei , ei = \ + \j\ + 2a2n2 ■ ( 3 - 3 ) 
At large n, corresponding to large Q, this implies 

a{z) ~ k lZ ^ an . (3.4) 
The function h can then be found from the small z limit of Eq. ( I2.59P 

h( z ) ~ ^ = V^ahz^- 1 . (3.5) 
n 

3.2 Large z asymptotic behavior 

Far from the monopole, a finite energy solution requires that a tends to n. It will prove 
convenient to define the function 

j3(z) = n — a(z) , (3.6) 
which tends to at large z. Eq. (12.6 ip may be re-expressed in terms of (3 as 

We will now assume that may be expanded at large z such that the leading term is 

P( z ) = k 2 z e2 e- mz (3.8) 

where m > 0. Substituting this into (13. 7p one finds 

eo — a 2 n 2 

+ ■ ■ ■ = , (3.9) 



,2 
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where the ellipsis denote exponentially suppressed terms, which may be canceled by sub- 
dominant corrections to (13. 8p . Therefore we demand only that the first term vanishes, 
fixing C2 = a 2 n 2 . 

Again the Higgs field profile function h is easily found from (I2.59P 

BPS monopoles are characterized by the boundary condition that h tends asymptotically 
to the vacuum expectation value gv (recall that we have rescaled the coupling into the 
scalar field). Therefore m is determined entirely by this boundary condition 

m = gv . (3-11) 

Physically, m is the mass of the W^-bosons, and the exponential decay in Eq. (13.81) is just 
that of a massive field. Our final asymptotic form for the function /3(z) is hence 

p(z) ~ k 2 z a%n2 e- 9VZ . (3.12) 



3.3 Connecting the regimes and the monopole bag 

We have studied two regimes: one, at small z, where a <C n and another, at large z, where 
(3 = n — a <C n. Clearly it is important to determine at what value z = z the solution of 
(I2.6ip interpolates between these two regimes. Neither approximation of z leads to a well- 
controlled expansion in the intermediate regime. However one may arrive at a qualitative 
understanding of the solution by imagining that both expansions are roughly correct, in 
the sense that will be described below, at Zq. This is potentially a dangerous assumption, 
but numerically we have verified that the results of this subsection are indeed correct. 

We will define z to be the midpoint of the function a 

Tl 

a(z ) = 0(z o ) = - . (3.13) 



Therefore the matching of the two limits yields 

a'(z ) (3'(z ) 



(3.14) 



a(z ) /3(z ) 

This relation is exact. But the consequences may be approximated by substituting the 
asymptotic behaviors (13. 4p of a at small z and (13.121) of ft at large z which gives an 
estimate for 

an f r- \ a 2 n 2 

zq [V2 + an)~ , 3.15 

gv \ J gv 
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where in the rightmost expression we have approximated n 1. From this expression it 
is clear that the 1/z tail of the Higgs field gives the monopole its size oc n 2 . 

The value of z in Eq. (I3.15P is the approximate size of the monopole, the value of the 
radius at which the H^-bosons begin to fall exponentially. At charge 1, corresponding to 
n — 1, one recovers the fact that the monopole size is the inverse W^-boson mass. However 
more generally it reveals, as was conjectured in Ref. [5], that the radius of the monopole 
is proportional to the charge Q = n 2 . 

One may also determine the thickness w of the boundary region, the distance over 
which a tends from a value near to a value near n. This depends on how quickly a 
changes at the boundary z 

a'(zo) = V2an^- ~ (3.16) 

valid for n ^> 1. The distance over which a changes by n units is then approximately 

n V2an 

w ~ — - — r ~ . (d.17) 

a'(zo) gv 



Therefore the width of the boundary of the monopole is proportional to n = y/Q, while 
the radius of the monopole is proportional to n 2 = Q. Thus at large n the walls of 
the monopole are much thinner than its radius, confirming a conjectured description of 
these solutions as monopole bags in Ref. [5]. However one should note that, given the 
large moduli space of solutions, it seems possible that there are other sequences of BPS 
solutions with approximate spherical symmetry in the large Q limit but with different 
qualitative radial profiles, and so the bag description may not apply to them. 

The boundary of the monopole is not only relatively thin, with a width w much smaller 
than z at large n, but also it is very sharp, as exponentially decays at large z as seen 
from Eq. (13. 12p . As we have mentioned, this is due to the fact that the gauge field is 
massive. On the other hand, the massless Higgs field decays only as 0(l/z), as seen 
in Eq. (I3.10p . Therefore the Higgs field does not exhibit such a sharp transition at the 
boundary of the monopole. 



4 Numerical results 

In this section we will provide a few numerical solutions to the ODE (I2.6ip illustrating 
the radial profile functions for the monopole. We have used a shooting method to find the 
solution using Eq. (13. 4 p as the initial condition and hence ki as the shooting parameter. 
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A solution with Q = 81, n — gv — 9, a — 1.9 is shown in Fig. [T] and its energy density 
is shown in Fig. [2J This solution has a total energy which is 54% higher than the BPS 
bound 4nvn 2 /g. 




10 20 30 40 50 100 200 300 400 

z z 



Figure 1: Left panel: the profile function a. Right panel: the profile function h both as 
functions of z for Q — 81, n = gv — 9, a — 1.9, zq = 25.8 and a = n/a. 



Energy density 




z z Tl 

Figure 2: Left panel: the radial energy density f° drj rjg 2, H cone as function of z for the 
monopole with Q = 81, n = vg = 9, a = 1.9. Right panel: the energy density g 2 H cone of 
the cone. 

In order to check the size of the monopole estimated in Eq. (I3.15p . we have numerically 
calculated the value of Zq such that a(z ) = fi/2 as shown in Fig. El 

Finally we have checked the excess energy of the solutions compared to the BPS- 
saturated ones which is shown in Fig. SJ The solutions turn out to exceed the BPS bound 
by roughly 54% at large n. 
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Figure 3: The approximate monopole size gvzo which is predicted to be In (y/2 + In) 
and shown as the red line (we have set o = 1.9). The numerical results are shown with 
their corresponding errorbars. 
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Figure 4: The total excess energy of the monopole compared to the BPS-saturated energy 
percentage as a function of n — y/Q. It is seen that the energy excess is roughly constant 
at large n around 54%. 
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5 Future directions 

First of all we should provide a word of caution due to the fact that the BPS solutions 
do not share the factorization property of our configurations. For this reason we have 
calculated the radial profiles in terms of angular expansion parameters expanded around 
the origin. This means that our monopole configurations are only approximately BPS, 
and in particular do not provide time-independent solutions of the equations of motion. 
However we do believe they capture some key features of those true BPS monopoles 
which are spherically symmetric. As a support to this claim we found numerically that 
at high Q the total energies of our configurations exceed that of true BPS solutions by a 
Q-independent factor of roughly 54%. 

As we have mentioned, it may seem as though we are studying a particularly diffi- 
cult and uninteresting part of the moduli space of solutions. We would however like to 
conjecture that in a certain class of theories it is the most interesting part: 

Conjecture 1 The approximately spherically symmetric BPS monopoles are the only 
ones which survive the strong non-BPS deformation described below. They all reduce to 
the same non-BPS configuration. 

Our deep interest lies in non-BPS monopoles in which a Higgs potential is included 
for the scalar We are interested in these monopoles because of a series of perhaps 
coincidental facts relating Q — 1 non-BPS monopoles with the dark matter halos of the 
many minimal size dwarf spheroidal galaxies which have recently been discovered in our 
local group, for example by the Sloan Digital Sky Survey. Some of the most striking 
similarities, which in fact are shared by other dark matter dominated galaxies such as 
larger dwarf and low surface brightness spiral galaxies, are as follows: 

1) Dark matter halos, like topological solitons, have a minimum mass. For solitons this 
corresponds to the charge Q — 1. The lightest satellites of the Milky Way have masses 
of about 1O 7 M within 1,000 light years of their center [TTJ and between 10 7 and 10 8 M 
within 2,000 light years [12J. This observed minimum dark matter mass leads to the dwarf 
galaxy problem: For particle models of dark matter, like WIMPs, consistency with this 
minimum mass is a problem because simulations generally suggest 4 to 400 times more 
dwarf galaxies in our local group than have been observed, most of which should be much 
lighter than the minimum observed mass. In fact, only one dwarf galaxy (ComaBer [13] ) 
has been observed with a mass near 1O 6 M , and it is very elongated and irregular and 
appears to be in the process of being ripped apart by tidal forces. 
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The existence of a topological charge, equal to one, for these small galaxies not only 
explains the fact that no small dwarf galaxies have been seen, but also the related fact 
that smaller dark matter bodies cannot exist. In this way one avoids the fatal gravitational 
lensing constraints faced by other MACHO dark matter models. Indeed the upper limit 
of the range of radii of dark matter candidates excluded by gravitational lensing is many 
orders of magnitude smaller than these 1000 light year solutions, and so these monopoles 
are too large to be excluded by the lensing bounds. 

2) At least in cases in which there is enough visible matter to determine the density profile, 
dark matter dominated galaxies and non-BPS monopoles have cores with relatively con- 
stant densities, with intermediate regions with 1/r 2 densities and external regions with a 
faster radial fall-off. For Q — 1 non-BPS monopoles this 1/r 2 intermediate region density 
profile is inevitable, unlike the the higher Q BPS 1/r 2 density profile in the Bolognesi 
galaxy bags of Ref. [S] in which there are many inequivalent choices of r-dependence, such 
as solutions which the author called planets, etc. 

3) The cores of these non-BPS solutions in many cases naturally contain black holes p~H 
[T5| dnj [T7]. In the case of the galaxies, models often suggest that there has not been 
enough time to form the supermassive black holes known to inhabit most galactic cores, in 
addition there are even some claims of supermassive black holes without luminous galactic 
hosts. These problems are both naturally explained if the black hole is an integral part 
of the monopole solution, as it is in many models. The gradual consumption of stars, 
gas and dark matter particles is, in this scenario, no longer the main mechanism driving 
supermassive black hole formation. 

4) The simplest model in which non-BPS 't Hooft-Polyakov monopoles exist is a Georgi- 
Glashow modefl with a simple Abelian Higgs quartic potential. In this case, a Q = 1 
non-Abelian monopole with the radius r and mass M of the smallest dwarf galaxies arises 
if the value of the Higgs VEV is about v ~ ^ho>Mjr ~ 10 14 GeV. This number only 
changes by about a factor of 2 depending on whether the luminous region is within the core 
or the intermediate radius regime. Had v been above the Planck scale, gravity would have 
dominated over the Georgi-Glashow interactions and the whole solution would have been 
a hairy black hole instead of a dwarf galaxy. Had it been smaller than about 100 eV, these 
monopoles could not have formed in time for dark matter to have played its crucial role in 
the oscillations of primordial plasma which reproduces the oscillation spectrum observed 

4 Here we are considering a new sector, these gauge symmetries have no obvious relation with standard 
model or GUT symmetries and we do not specify the charges of standard model particles under this new 
gauge symmetry. 



21 



in the CMB. Given that the two physical inputs in this calculation are of galactic scales, 
the fact that the output is a particle physics scale in this relatively narrow acceptable 
window is for us miraculous. If one naively uses the rotation curves of slightly larger 
dwarves one may similarly conclude that the Higgs coupling is of order A ~ 10~ 97 , had 
it entered v with a different power, even a fourth root, the relation with dwarf galaxies 
would have been ruined. 

5) Similarly the 1,000 light year scale radii of these solutions imply that they form when the 
universe is about 1,000 years old [18]. Again, this is in time to help increase the intensity 
of fluctuations in the primordial plasma as is required by observations of the CMB. Had 
dwarf galactic radii been larger by a factor of 100, they would have formed too late and 
an inconsistency would have arisen. 

6) While the monopole core excluding gravity has a constant density, and with gravity may 
host a black hole, the core itself is nonsingular. More precisely it avoids the cusp prob- 
lem of the ACDM model, in which many simulations predict galactic mass distribution 
profiles, such as the historic Navarro, Frenk and White profile [19], with density cusps in 
their cores, in stark contrast with observations. 

While these similarities are very strong, there is a serious problem with galaxy sized 
non-BPS Bolognesi bags as a dark matter candidate. Non-BPS monopoles repel, and 
as v is less than the Planck scale this repulsion would dominate over gravity and all 
galaxies would be minimal dwarves and would repel one another. There is a similar 
problem of course for visible matter, which is mostly made of protons which also repel. 
In the case of protons, the solution is quite complicated. First of all there are electrons 
which screen the interactions between protons. While electrons have antiparticles which 
have the same charge as protons, for reasons which have not yet been quantitatively 
explained by any model, there was a primordial excess of negatively charges electrons 
and positively charged protons. They did not annihilate each other because they carry 
different conserved charges. They can combine, forming hydrogen bound states or via 
inverse beta decay they can even merge beyond recognition into neutrons. However, due 
to the choice of parameters in the standard model, the later possibility is kinematically 
disfavored in the conditions that have existed in most of the universe since baryogenesis. 

We would like to propose that repulsion between non-BPS monopoles is avoided in a 
similar manner. Additional conserved charges are easily introduced in a Georgi-Glashow 
model by including charged fermions, which via the Jackiw-Rebbi mechanism provide an 
additional charge for each kind of monopole. If one adds two species of fermions, then 
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there are two kinds of charge, which can play a role analogous to baryon number and 
lepton number. Monopoles of different charges can have very different masses, in fact in 
M = 1 supersymmetric models some flavors are usually massless while some are massive. 
One can then demand that the dark matter halos are made of very massive magnetically 
positively charged monopoles which carry one kind of flavor charge, and that the screening 
is caused by light negative monopoles with the other flavor charge. This eliminates the 
problem of galaxies repelling another. 

But one still needs to worry about the stability of Q > 1 monopoles. These will 
be held together by gravity. Due to the mass of the scalar field <£>, the repulsive SU(2) 
dynamics will dominate over the attractive scalar dynamics at large distances, leading 
to a net repulsion. The gravitational interaction in general is insufficient to counter this 
repulsion, as v is much less than the Planck mass. However at large distances one expects 
the screening to play a role. Unfortunately the exact role played by this screening is 
highly model dependent and is not clear whether there exist any models which screen 
this repulsion sufficiently to allow it to be dominated by gravitational attraction, without 
adding any new interactions (playing the role of strong interactions in the proton analogy 
described above). Strong constraints on models also arise from the fact that one does not 
want the light monopoles to combine with the heavy monopoles, analogously to inverse 
beta decay, as the resulting bound state may not share the attractive features of the 
massive monopoles described above. 

Therefore the model-independent predictions for Q > 1 monopoles are limited by 
ambiguities in the screening mechanism. Nonetheless a number of very firm predictions 
can be made. First of all, just as the dwarf galaxy problem is a gap between the mass of 
globular clusters at Q = and dwarf galaxies at Q — 1, there must also be a gap between 
Q = 1 and Q = 2. This prediction is much more general than the monopole dark matter 
proposal discussed in this section, but extends to any topological soliton dark matter 
candidate which solves the dwarf galaxy problem by identifying minimal spherical dwarf 
galaxies with Q = 1 solitons. The masses of these galaxies are at best known at the 100 
percent level, and so with current data such a gap cannot be verified. However one may 
hope that radio surveys of gas in our galactic neighborhood such as that which will be 
performed by the FAST telescope starting in 5 years will be able to test this claim. 

Another model independent prediction is that, while the charge Q is determined by 
the flat part of the galactic rotation curve, the radius of the core must be proportional 
to and, even more surprisingly, the outer radius of the region with the flat rotation 
curves must at large Q be nearly Q-independent. This counter- intuitive prediction may 
already rule out these models, as it requires spatial extents of dwarf galaxy dark matter 
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halos to extend far beyond their most distant stars, but it is necessary for the convexity of 
the galactic mass as a function of Q, which in turn is necessary to prevent these galaxies 
from exploding. 

Acknowledgments 

We are eternally grateful to Stefano Bolognesi and Malcolm Fairbairn for many useful and 
enlightening discussions. JE is supported by the Chinese Academy of Sciences Fellowship 
for Young International Scientists grant number 2010Y2JA01. SBG is supported by the 
Golda Meir Foundation Fund. 

References 

[1] R. S. Ward, "A Yang-Mills Higgs Monopole of Charge 2," Commun. Math. Phys. 
79, 317-325 (1981). 

[2] M. K. Prasad, P. Rossi, "Construction Of Exact Yang-mills Higgs Multi - Monopoles 
Of Arbitrary Charge," Phys. Rev. Lett. 46, 806 (1981). 

[3] M. K. Prasad, P. Rossi, "Rigorous Construction Of Exact Multi - Monopole Solu- 
tions," Phys. Rev. D24, 2182 (1981). 

[4] M. K. Prasad, "Exact Yang-mills Higgs Monopole Solutions Of Arbitrary Topological 
Charge," Commun. Math. Phys. 80, 137 (1981). 

[5] S. Bolognesi, "Multi-monopoles, magnetic bags, bions and the monopole cosmological 
problem," Nucl. Phys. B752, 93-123 (2006). |hep-th/0512133] . 

[6] K. -M. Lee, E. J. Weinberg, "BPS Magnetic Monopole Bags," Phys. Rev. D79, 
025013 (2009). |arXiv:0810.4962l [hep-th]]. 

[7] D. Harland, "The Large N limit of the Nahm transform," jarXiv: 1102.30481 [hep-th]]. 

[8] N. S. Manton, "Monopole Planets and Galaxies," jarXiv: 1111 .29341 [hep-th] ] . 

[9] Stefano Bolognesi, Private communication in 2010. 

[10] P. Sutcliffe, "Monopoles in AdS," JHEP 1108 (2011) 032. jarXiv: 1 104. 18881 [hep-th] ] . 



24 



[11] L. E. Strigari, J. S. Bullock, M. Kaplinghat, J. D. Simon, M. Geha, B. Willman 
and M. G. Walker, "A common mass scale for satellite galaxies of the Milky Way," 
Nature 454 (2008) 1096 [arXiv:0808.3772l [astro-ph]]. 

[12] L. E. Strigari, J. S. Bullock, M. Kaplinghat, J. Diemand, M. Kuhlen and 
P. Madau, "Redefining the Missing Satellites Problem," Astrophys. J. 669 (2007) 
676 jarXiv:0704.18T7l [astro-ph]]. 

[13] V. Belokurov et al. [SDSS Collaboration], "Cats and Dogs, Hair and A Hero: 
A Quintet of New Milky Way Companions," Astrophys. J. 654 (2007) 897 
|arXiv:astro-ph/0608448] . 

[14] K. M. Lee, V. P. Nair and E. J. Weinberg, "A Classical instability of Reissner- 
Nordstrom solutions and the fate of magnetically charged black holes," Phys. 



Rev. Lett. 68 (1992) 1100 arXiv:hep-th/9111045|; K. M. Lee, V. P. Nair and 



E. J. Weinberg, "Black holes in magnetic monopoles," Phys. Rev. D 45 (1992) 2751 
jarXiv : hep-th / 9 112 008 1 ; A. Lue and E. J. Weinberg, "Magnetic monopoles near the 
black hole threshold," Phys. Rev. D 60 (1999) 084025 [arXiv:hep-th/9905223] . 

[15] P. Breitenlohner, P. Forgacs and D. Maison, "Gravitating monopole solutions," Nucl. 
Phys. B 383 (1992) 357; P. Breitenlohner, P. Forgacs and D. Maison, "Gravitating 
monopole solutions. 2," Nucl. Phys. B 442 (1995) 126 [arXiv:gr-qc/9412039| . 

[16] B. Hartmann, B. Kleihaus and J. Kunz, "Gravitationally bound monopoles," Phys. 
Rev. Lett. 86 (2001) 1422 [arXiv:hep-th/0009195] ; B. Hartmann, B. Kleihaus and 
J. Kunz, "Axially symmetric monopoles and black holes in Einstein- Yang- Mills-Higgs 
theory," Phys. Rev. D 65 (2002) 024027 |arXiv:hep-th/0108129| . 

[17] S. Bolognesi, "Magnetic Bags and Black Holes," jarXiv: 1005.46421 [hep-th]]. 

[18] T. W. B. Kibble, "Topology Of Cosmic Domains And Strings," J. Phys. A 9 (1976) 
1387. 

[19] J. F. Navarro, C. S. Frenk and S. D. M. White, "The Structure of Cold Dark Matter 
Halos," Astrophys. J. 462 (1996) 563 |arXiv:astro-ph/9508025] . 



25 



